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Abstract 

We perform a theoretical study based on dispersion relations of the reaction 77 — > 7r°7r° 

emphasizing the low energy region. We discuss how the /o(980) signal emerges in 77 — * tttt 

within the dispersive approach and how this fixes to a large extent the phase of the isoscalar 

S-wave 77 — * tttt amplitude above the KK threshold. This allows us to make sharper 

predictions for the cross section at lower energies and our results could then be used to 

distinguish between different tttt isoscalar S-wave parameterizations with the advent of new 

precise data on 77 — ► 7r°7r°. We compare our dispersive approach with an updated calculation 

employing Unitary Chiral Perturbation Theory (UxPT). We also pay special attention to 

the role played by the a resonance in 77 — > -kit and calculate its coupling and width to 77, 

O for which we obtain T(a -> 77) = (1.68 ± 0.15) KeV. 
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1 Introduction 



The study of the reaction 77 — > 7r°7r° offers the possibility of having two neutral pions in the 
final state as a two body hadronic system. In this way, the interactions among pions can be 
unraveled more cleanly than in other processes where three or more hadrons appear in the final 
state. In this respect, the K®2 decays also share this property but limited to energies below the 
kaon mass. Furthermore, the two photon fusion also offers the interesting property of coupling 
mostly to the charged components of the hadronic states, while strong interactions are blind to 
that. Therefore, the study of 77 fusion into two hadrons displays new information on the inner 
constituents, beyond that obtained in direct strong scattering data. In addition, 77 — > 7r°7r° has no 
either Born term, shown in fig. 1 for charged mesons, so that final state interactions dominate this 
reaction. In this respect, one can think of the interest of precise data on this reaction to distinguish 
between present low energy parameterizations of the nn isospin (I) S-wave, to further study the 
nature of the a resonance as well as that of the /o(980), and having another way to constraint 
their pole positions. The presence of the a and its tight relation to chiral symmetry was clearly 
established from those non-perturbative methods [1-5] that resum unitarity in the series of Chiral 
Perturbation Theory (%PT) [6, 7]. These methods were the first to show that the presence of a a 
pole is not incompatible with the non-linear a models. They also showed a clear convergence in 
the a pole position when using either lowest order or one or two loop ;\PT in the resummation, 
this is discussed in detail in ref.[5]. Recently, its associated pole was determined by solving Roy 
equations together with constraints from %PT to two loops [8]. Other recent determinations are 
refs.[9,10]. 

The reaction 77 — > 7r°7r° is also an interesting ground test for %PT [6, 7], since at lowest order 
is zero and at next-to-leading order (one loop) it is a prediction free of any counterterm [11, 12]. 
The one loop %PT calculations were performed before the advent of data but the prediction just 
above threshold departs rapidly from experiment [13] and only the order of magnitude was rightly 
foreseen. A two loop %PT calculation was done in refs.[14, 15], where the agreement with data 
[13] at low energies was improved. The three counterterms that appeared at 0(p 6 ) were fixed 
by the resonance saturation hypothesis [16]. Other approaches supplying higher orders to one 
loop xPT by taking into account unitarity [17] plus also crossed exchange resonances [18, 19] 
followed. Ref. [19] is a dynamical approach, with higher predictive power than dispersion relations 
(which make use of extensive experimental input from other related processes), though more model 
dependent. Ref. [19] is a Unitary (UxPT) calculation in production processes [19-24] and 

was able to simultaneously describe 77 — > 7r°7r°, 7r + 7r~, rjn , K + K~ and K°K° from threshold 
up to rather high energies, s 1 / 2 < 1.5 GeV. This approach was also used in ref. [25] to study the 
r] — > 7T°77 decay. 

Recently, ref. [26], making use of the approach of refs.[27, 28], focused on the low energy part 
of the reaction 77 — > ir°n°. Ref. [26] emphasizes the calculation of the 77 coupling of the a as 
a way to deepen in the understanding of its nature. This approach was extended in ref. [29] so 
as to largely remove the sensitivity for lower energies, ^/s < 0.8 GeV, on the uncertainty in the 
not precisely known phase of the 1 = S-wave 77 — > nn amplitude above the KK threshold. 
The novelty was to include a further subtraction in the dispersion relation for the 1 = S-wave 
77 — > 7T7T, together with an extra constraint to fix the additional subtraction constant. This was 
motivated by the use of an improved 1 = S-wave Omnes function in the lines of ref. [30]. We 
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shall here extend the findings of ref.[29] and discuss in detail how the / (980) peak, clearly seen 
recently in 77 — > tx + tx~ [31], can be generated within the dispersive method of ref.[29]. We also 
show how this can be used to fix rather accurately the phase of the 1 = S-wave 77 — > txtx 
amplitude above the KK threshold. As a result, the remaining ambiguity in this phase by tx in 
ref. [29] is removed. Interestingly, this allows to sharpen the prediction of the 77 — ■> 7r°7r° cross 
section for y/s < 0.8 GeV, up to the onset of D-waves, as compared with ref. [29]. This could then 
be used to constraint further different parameterizations of the low energy txtx I = S-wave. We 
also offer explicit expressions for several axial-vector and vector crossed exchange contributions 
to 77 — > Tin, that are used in the present analysis and that of ref. [29]. The Unitary %PT (UxPT) 
calculation of ref. [19] on 77 — > txtx will be reviewed and extended, so as to incorporate the same 
vector and axial-vector exchange contributions as in the dispersive analysis. We finally discuss 
the a coupling to 77 and determine Via — > 77) = (1.68 ±0.15) KeV. Other papers devoted to the 
calculation of the two photon decay widths of hadronic resonances are [32-43]. 

The content of the paper is as follows. In section 2 we discuss the dispersive method of ref. [27] 
and its extension by ref. [29] for calculating 77 — > txtx. Here we pay attention to the intrinsic 
difficulties of the approach of refs. [26-28] to generate the / (980) peak in 77 — > txtx and how this 
is overcome by the approach of ref. [29] . In this way, one can use this extra information to improve 
the calculation at lower energies. In section 3 we discuss and revise the UxPT approach of ref. [19]. 
We also give in this section several formulae corresponding to the Born terms and exchanges of 
vector and axial-vector resonances to be used as input in both approaches. The resulting cross 
section (7(77 — > 7r°7r°) is discussed in section 4. The 77 decay width of the a resonance and its 
coupling to 77 is considered in section 4.3. Conclusions are elaborated in section 5. 

2 Dispersive approach to 77 — » tttt 

We first review the approach of refs. [26-28] to establish a dispersion relation to calculate the two 
photon amplitude to txtx. This approach has a large uncertainty above yfs > 0.5 GeV, which 
dramatically rises with energy such that for y/s = 0.5, 0.55, 0.6 and 0.65 GeV it is 20, 45, 92 
and 200%, see fig.3 of ref. [26]. This uncertainty was largely reduced by the extended method 
of ref. [29]. We discuss here how this method allows to disentangle also the role of the / (980) 
resonance in 77 — > txtx and the way that this can be used to sharpen further the results of ref. [29] . 

2.1 Dispersion relation and Omnes function 




Figure 1: Born term contribution to 7(^1)7(^2) — ► P + (pi)P (P2). 
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Due to the null charge of the tt° there is no Born term, shown in fig.l, for 77 — > 7r°7r°. As 
remarked in ref . [28] , one then expects that only the S-wave would be the important partial wave 
at low energies, y/s < 0.6 — 0.7 GeV. For 77 — > ir + ir~, where there is a Born term because of 
the exchange of charged pions, the D-waves have a relevant contribution already at surprisingly 
low energies due to the smallness of the pion mass. In the following, we shall restrict ourselves 
to the S-wave contribution to 77 — > tttt, being particularly interested in the low energy region. 
Explicit calculations [19] show that the D-wave contribution in the 77 — > 7r°7r° cross section at 
y/s ~ 0.65 GeV amounts less than a 10%, and this rapidly decreases for lower energies. Let us 
consider the S-wave amplitude 77 — > (7T7r)j, Fj(s), where the two pions have definite / = or 
2. Notice that charge conjugation of the photons is even and this, together with Bose symmetry, 
excludes odd orbital angular momentum (and hence odd isospins) for the two pions. The function 
Fj(s) on the complex s— plane is analytic except for two cuts, the unitarity one for s > 4m 2 and 
the left hand cut for s < 0, with m n the pion mass. Let us denote by £/(s) the complete left hand 
cut contribution of Fj(s). Then the function F I (s) — Lj(s), by construction, has only right hand 
cut. 

Refs. [26-28] consider the Omnes function u>i(s), 



oui(s) = exp 



TT J 4m 2 S'(s' - S 



(2.1) 



with (f>i(s) the phase of Fj(s) modulo tt, chosen in such a way that (f>i(s) is continuous and 
0/(4m^.) = 0. This function was criticized in ref. [30]. The 1 = S-wave phase shift, 5 n (s) , has 
a rapid increase by tt around sk = ^m^, with m K the kaon mass, due to the narrowness of the 
/o(980) resonance on top of the KK threshold. Let us denote by ip(s) the phase of the tttt — > tttt 
1 = S-wave strong amplitude such that it is continuous, modulo tt when crossing a zero, and 
(p(4m%) = 0. This phase is shown in fig. 2 together with 5 n (s)o and S n (s)2, where the latter is the 
tttt I = 2 phase shift. Now, if one uses <p(s) instead of o (s) in eq.(2.1) for illustration, the function 
ujq(s) is discontinuous in the transition from S n (s K ) — > tt — e to S n (sK)o — > n + e, e — > + . In the 
former case |^o(s)| has a zero, while in the latter it becomes +00. This discontinuity is illustrated 
in fig. 3 by considering the difference between the dashed and dot-dashed lines, respectively. This 
discontinuous behaviour of ujq{s) under small (even tiny) changes of 5 7r (s)o around sk was the 
reason for the controversy regarding the value of the pion scalar radius (r 2 )J between refs. [44-46] 
and ref. [47]. This controversy was finally solved in ref. [30] where it is shown that Yndurain's 
method is compatible with the solutions obtained by solving the Muskhelishvili-Omnes equations 
for the scalar form factor [48-50]. The problem arose because refs. [44, 45] overlooked the proper 
solution and stuck to an unstable one. 

Given the definition of the phase function <f>i(s) in eq.(2.1) the function Fj(s)/uj(s) has no 
right hand cut. Then, refs. [26-28] perform a twice subtracted dispersion relation of (Fj(s) — 
L i (s))/uj i (s), 

Fj(s) = Lj(s) + aj uj(s) + cjs ujj(s) + -uj (s) / J ) M ds " (2 ' 2) 

TT J Am 2 S' {S' - S)|W/(S')| 

Low's theorem [51] implies that Fj — > Bj(s) for s — > 0, with Bj the Born term contribution, shown 
in fig.l. If we write Lj = Bj + with Rj — * for s — > 0, as can always be done, then the Low's 
theorem implies also that Fj — L/ — > for s — > and hence a/ = 0. 
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Figure 2: The phase shifts 5,r(s)o and ^(5)2 and the phase (p(s). Experimental data are from refs. [52-56] 
for 1 = and refs. [57, 58] for 1 = 2. The insert is the comparison of CGL [59], dot-dashed line, PY [60], 
dashed line, and AO [10], double-dash-dotted line, with the accurate data from K e ^ [55, 56]. For AO we 
have not shown its error bar in order not to overload the figure, see fig.l of ref.[10] where it is given. 

For the exotic 1 = 2 S-wave one can invoke Watson's final state theorem* 1 so that 02 (s) = 
5 7r (s) 2 . For 1 = the same theorem guarantees that 4>o(s) = 5 n (s)o for s < 4m K . Here one 
neglects the inelasticity due to the Ait and 6tc states below the two kaon threshold, an accurate 
assumption as indicated by experiment [52,53]. Above the two kaon threshold sk = 4m^, the 
phase function <po{s) cannot be fixed a priori due the onset of inelasticity. This is why ref.[26] 
took for s > sk that either 4>o(s) ~ 5 n (s)o or (j)o(s) ~ 5tt(s)o — tt, in order to study the size of 
the uncertainty induced for lower energies. It results, as mentioned above, that this uncertainty 
becomes huge already for > 0.55 GeV [26]. 

Let us now explain how this problem was settled in ref.[29]. Inelasticity is again small for 
1-1 ~ y/s ^ 1-5 GeV [52,53]. As remarked in refs. [30, 44], one can then apply approximately 
Watson's final state theorem and for Fq(s) this implies that (fio(s) ~ modulo tt. Here 

5( + )(s) is the eigenphase of the nn, KK 1 = S-wave S-matrix such that it is continuous and 
8^ + \sk) = 5 n (sK)o- in refs. [30, 45] it is shown that <5( + )(s) ~ 5,r(s)o or 5 7r (s) — it, depending 
on whether 5 7r (sx)o > ^ or < tt, respectively. In order to fix the integer factor in front of n in 
the relation (fio(s) ~ 5^~\s) modulo tt, one needs to devise an argument to follow the possible 
trajectories of 4>o(s) m the narrow region 1 < \/s < 1.1 GeV, where inelasticity is not negligible. 

^ 1 This theorem implies that the phase of F](s) when there is no inelasticity is the same, modulo tt, as the one 
of the isospin / S-wave tttt elastic strong amplitude. 
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Figure 3: |wo(s)| from eq.(2.1) with 5- k (sk)o < dashed-line, and S- k (sk)o > 7r, dot-dashed line. The 
solid line corresponds to the use of ^Iq(s) in eq.(2.3) for the latter case. Here ip(s) is used as 0o(s) in 
eq.(2.1) for illustrative purposes. 

The remarkable physical effects happening there are the appearance of the /o(980) resonance on 
top of the KK threshold and the cusp effect of the latter that induces a discontinuity at sk in 
the derivative of observables, this is clearly visible in fig. 2. Between 1.05 to 1.1 GeV there are 
no further narrow structures and observables evolve smoothly. Approximately half of the region 
between 0.95 and 1.05 GeV is elastic and (f)o(s) = 6^(3)0 (Watson's theorem), so that it raises 
rapidly. Above 2rriK — 1 GeV and up to 1.05 GeV the function <po(s) can keep increasing with 
energy, like 5 n (s)o or (/?(s) for S n (sK)o > tt, and this is always the case for the corresponding 
phase function of the strange scalar form factor of the pion [30]. It is also the behaviour for <fio{s) 
corresponding to the explicit calculation of ref.[19]. The other possibility is a change of sign in the 
slope at sk due to the KK cusp effect such that 0o( s ) starts a rapid decrease in energy, like ip(s) 
for ^(5^)0 < 7r, fig.2. Above ^/s = 1.05 GeV, 0o( s ) matches smoothly with the behaviour for 
y/s > 1.1 GeV, which it is worth recalling that it is constraint by Watson's final state theorem. As 
a result of this matching, for ^/s > 1 GeV either (f>o(s) ~ S n (s) or0 o (s) ~ 5 n (s) — it, corresponding 
to an increasing or decreasing 4>q(s) above sk, in order. Our argument also justifies the similar 
choice of phases in ref.[26] above sk to estimate uncertainties. 

Let us define the switch z to characterize the behaviour of <po(s) for s > sk, and close to sk, 
such that z — +1 if <fro(s) rises with energy and z — —1 if it decreases. Let Si be the value of s at 
which 0o(si) = 7T- As shown in ref.[30], instead of eq.(2.1), one should better consider the Omnes 
function, 



with 9(z) — 1 for z — +1 and for z — — 1. Because of the first order polynomial in fio(s); with 
a zero at s±, one can guarantee a continuous behaviour for Qq(s) when moving from z — — 1 to 
z — +1, passing from a deep minimum to a zero. This is shown in fig.3 by comparing the solid 
line with the dashed one, both computed with eq.(2.3) and using ip(s) as 0o( s ) f° r illustrative 
purposes. Both lines are very close to each other as corresponds to a continuous transition. Note 
that Si is the only point around sk where the imaginary part of u)$(s) vanishes, and this fixes 
unambiguously the position of the zero in eq.(2.3), as discussed in ref.[30]. For 1 = 2 since the 
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phase shifts are small and smooth [57,58], see fig. 2, the issue of the discontinuity under changes 
in the parameterizations does not arise and we keep using 0*2 (s), eq.(2.1). 

Now, we perform a twice subtracted dispersion relation for (F (s) — L (s))/^ (s), employing 
ilo(s) instead of c<j (s). Instead of eq.(2.2) one obtains 

f (s) = l (s) + c oS n (s) + -no (a ) s , 2{s ,_ s)msl)l ds + ew^^Fofr) - L ( Sl )) . 

(2.4) 

In the previous equation we introduce <f>o(s) that is defined as the phase of ilo(s). Let us note 
that in the case z = +1 the phase of fi (s) f° r s > Si is not 0o( s ) but <po(s) — n, due to the factor 
1 — (s + ie)/si in fio(s), ec i-(2-3), present in this case. For / = 2 we use eq.(2.2). 

It is worth mentioning that eq.(2.4) for / = and z = +1 is equivalent to perform a three times 
subtracted dispersion relation for (F (s) — L (s)) / cu (s) , instead of the twice subtracted dispersion 
relation in eq.(2.2) as proposed in refs. [26-28]. In eq.(2.4), two subtractions are taken at s = 
and a third one at si. We could have taken the third subtraction at s = as well, so that for 
z — +1 one has 

F s = L (s) + c sw (s) + d s 2 w (s) + ^ / °\ ' J°) ds' , (2.5) 

J Am 2 S'-\S' - S)\U) Q (S')\ 

Nevertheless, we find more convenient eq.(2.4) as it is physically motivated by the use of the Omnes 
function eq.(2.3) that is continuous under changes in the parameterization of the 1 = S-wave 
S-matrix. Furthermore, eq.(2.4) is more suited to constraint the /o(980) energy region in order to 
implement the condition 3 below, because s± ~ 1 GeV 2 and it coincides very approximately with 
the mass of this resonance. 

Let us denote by Fn(s) the S-wave 77 — > 7i°n° amplitude and by Fc(s) the 77 — > 7r + 7r~ one. 
The relation between F , F 2 and F N (s), F c (s) is 

Fc(s) = ~F -^F 2 . (2.6) 

We are still left with the unknown constants cq, c 2 for / = and 2, respectively, and Fq(si) — 
L Q (si) for / = and z = +1. In order to determine them we impose the following conditions: 

1. F c (s) — B c {s) vanishes linearly in s for s — > and we match the coefficient to the one loop 
XPT result [11,12]. 

2. Fjv(s) vanishes linearly for s — > as well and the coefficient can be obtained again from one 
loop X PT [11,12]. 

3. For 1 = and z — +1 one has in addition the constant Fo(si) — L (si). Its value can be 
restricted taking into account that F N (s) has an Adler zero, due to chiral symmetry. This zero is 
located at sa = rn 2 in one loop %PT and moves to sa = 1.175m 2 in two loop %PT [14]. A rather 
large correction around 18% results, which prevents us from taking sa straightforwardly as given 
by xPT- In turn, we obtain that the value of the resulting cross section (7(77 — > 7r°7r°) around the 
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/o(980) resonance is quite sensitive to the value of F (si) — L (si). This is so because this constant 
appears in the last term in eq.(2.4) which dominates the energy region around the /o(980) since 
^o(si) = 0. Though the dispersive method is devised at its best for lower energies, it is also clear 
that it should give, at least, the proper order of magnitude for (7(77 — > 7r°7r°) in the /o(980) region. 
E.g., other models [19,61,62] having similar basic ingredients describe that energy region quite 
well indeed. We then restrict F (si) — £o( s i) so that (7(77 — > 7r°7r°) < 40 nb at si. Note that the 
D-wave component of this cross-section at si ~ 1 GeV, due to the tail of the / 2 (1270) resonance, 
is not longer small and hence the S-wave contribution should be significantly smaller than the 
experimental value < 40 nb. This is clearly shown by experiment [13], see fig. 9. The S- and 
D-wave contributions add incoherently in the total cross section because of the integration on the 
azimuthal angle as the /2(1270) resonance signal is overwhelming helicity 2 [27]. We have checked 
that this upper bound for the peak of the / (980) in 77 — > 7r°7r° is equivalent to impose that the 
77 width of the / (980), r(/ (980) -> 77), lies in the range 205tf 3 (stat)tn 7 7 (sys) eV, determined 
by the Belle Collaboration in ref.[31]. The width r(/o(980) — > 77) is calculated analogously as 
done in section 4.3 for the 77 width of the a resonance. We shall see that the effect of this rather 
large uncertainty allowed at 1 GeV, see fig. 9, is very mild at lower energies, as shown in fig. 6. 

It is worth discussing further the improvement achieved by eq.(2.4), first derived in ref. [29], in 
the understanding of the role played by the / (980) resonance in 77 — > -kit. The resulting |F (s)| 
given by eq.(2.2) from ref. [26] explodes for z — +1 around si, as |u>o(s)| becomes extremely large 
for this energy as shown in fig. 3 by the dot-dashed line. To avoid this large number the coefficient 
of Wq(si) in eq.(2.2) is required to almost vanish, 



Thus, a dramatic fine tuning of (po( s ) is implied, so as the integral above must almost cancel c . 
Note that Co is a real number and this integral is real only at s — s±, precisely where it is evaluated. 
Then the cancellation in eq.(2.7) is possible. But cq, following the approach of ref. [26, 28], is already 
determined by the requirement that the Adler zero for Fn(s) is at sa = m ^/2, rn% or 2m 2 [26]. 
The subsequent uncertainty in cq, together with the uncertainty in the knowledge of (f>o(s) above 
the KK threshold, invalidates in practice the fine tuning in eq.(2.7). This makes that the results of 
ref. [26] for (7(77 — > 7r°7r°) diverge rapidly with energy above y/s ~ 0.5 GeV for z = +1 as remarked 
above, see fig. 3 of ref. [26]. Refs.[18, 27, 28] did not discuss this source of large uncertainty in their 
calculations but, since these papers follow the same scheme as ref. [26], it is certainly present. One 
of the main results that follow from eq.(2.4) is to showing that by using VLq{s) one can cure such 
diverging results for z — +1. This issue is solved because for s around s\ one has that Qq(s) ~ 
and the diverging effects due to cj (s) for z = +1 are isolated in the last term and controlled by the 
factor Fq(si)—L (si). The small signal of the / (980) implicit in the 17(77 ~^ 7r°7r°) data of ref. [13], 
and explicitly seen in 77 — > 7t + tt~ [31], implies a vanishing value of the constant Fq(si) — Lq(si). 
In our actual calculations we obtain typically (F (si) — L (si)) / L (si) ^ 10~ 2 . One can indeed 
understand the smallness of this number by a continuity argument. For that, let us envisage a 
mathematical continuous transition from z = —1 to z = +1 by taking ip(s) as 0o( s ) i n eq.(2.4). 
The situation z = — 1 is accomplished when S w (sk)o < tt and z — +1 when ^(s^jo > n, see 
fig. 2. Now, imposing that F (s) is continuous under the transition S n (sK)o — > tt — e (where the 




(2.7) 
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last term in eq.(2.4) is absent) to ^(s^jo — > 7T + e, e — > + , where it is present, implies that 
Fo(si) — Lo(si) = for 6k{sk)o = 0. Since (f>o(s) is quite close to ip(s) both for z = +1 and — 1, 
the continuous transition considered is not just academic because S^sk) — n.* 2 As a result one 
should expect that F (s 1 ) — L (s 1 ) ~ 0. 

Since the /o(980) resonance gives rise to a small peak in the precise data on 77 — > 7t + tt~ [31], 
then 0o ( s ) must increase with energy above and the case with z = +1 is the one realized in 
nature. As a result 4>o(s) ~ ^(5)0, at least up to around < 1.5 GeV. Otherwise z = — 1 in 
eq.(2.4) and because |u;o(s)| nas a deep around the / (980) mass there is no a local maximum 
associated with this resonance in |-Fo( s )| but a minimum. One should keep in mind that z — +1 
and only this case will be considered in the subsequent. 

In ref. [29] both cases z — ±1 were discussed and kept, so that the differences in the results 
by taking one case or the other were included in the uncertainties of the approach, e.g. when 
presenting (7(77 — > 7r°7r°) for y/s < 0.8 GeV. In the latter reference the main focus was driven 
towards the drastic reduction of the uncertainty in the results for the two choices of (po(s) above 
sk- Here we are more accurate as we have singlet out the z — +1 case as the proper one for 4>o{ s ) 
and s > sk- 

We can estimate the uncertainty in the approximate relation 0o( s ) — <^r(s)o f° r 4m|- < y/s < 
1.5 GeV by similar arguments to those employed in ref. [30]. Let the irir and KK 1 = S-wave 
S— matrix, 

/ r)e 2iS " i^l-r] 2 e i( - s * +SK ) 

' ~~ V V 1 -r] 2 e i{ ^ +5K ) 7]e 2i5li 

with i] the elasticity parameter. This matrix is diagonalized with the orthogonal matrix C, 

C=( COs9 a S[n6 a ) . (2.9) 
y — sm cos J 

In eq.(2.8) 5^ = 5tt(s)o and 5k is the kaon phase shifts. Explicit expressions for cos# and sin^ 
in terms of the 1 = S-wave nn and KK phase shifts and the elasticity parameter are given in 
ref. [30]. In this reference one can also find the expression for the eigenvalues e 2t<5(+) and e 2%&( --*> . 
This diagonalization allows to disentangle two elastic scattering channels and then the two photon 
1 = S-wave amplitudes attached to every of these channels, Y\ and T' 2 , will satisfy the Watson's 
final state theorem in the whole energy range. Because of this diagonalization one has, 




F = costfr; + sin#r' 2 , 

F K = cos^-sin^r; . (2.10) 

Here Fk is the 1 = S-wave 77 — > KK amplitude. Now, when 77 — > 1 then sin# — > as 
e = a/(1 — rj)/2 [30] and F (s) is given approximately by r[(s). With tanOTyV^ = ee v , we 

#2 See ref. [45] for a detailed discussion on the value of S v (sk). 
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rewrite F above as, 



F (s) = cos6 T[ [1 + e(cosp + i sinp)] 

„ w /, x f , iesinp \ 

= cos#r 1 + ecosp 1 + — 

1 V 1 + ecosp J 

= cos6T[(l + ecosp)e iX + 0(e 2 ) . (2.11) 

with A = e sinp/(l + ecosp). From the previous equation it follows that 

Ms) = K(s)o + A + 0(e 2 ) (2.12) 

with e — > for rj — > 1. Using 77 = 0.8 in eq.(2.12) for the range 1.1 < y/s < 1.5 GeV (from the K- 
matrix in eq.(4.44) 77 ~ 1), one ends with |e| ~ 0.3. Taking into account that S n (s)o is larger than 
37r/2, fig. 2, then one has a relative correction of around a 10%. We can also proceed similarly in 
the region 1 < y/s < 1.1 GeV. Applying now eq.(2.12) with r\ > 0.6 from the nn — > KK scattering 
experiments [63], and taking S v (s)o > vr in the previous energy range, we end up with a ~ 20% 
uncertainty. We shall incorporate in our error analysis the previous uncertainties. 

We now consider the phases used above sh = 2.25 GeV 2 . For 1 = 2, where final state interaction 
give rise to small corrections, we directly take the extrapolation of our fits. Notice that due to 
the way that the vector and axial-vector exchanges are calculated from the chiral Lagrangians 
presented in the next section, one has that Lf(s) diverges linearly in s for s — > +00. On the other 
hand |a>/(s)| oc s~ Vl with i>j = 4>i{po)/n. Taking this into account for I = 2 the integrand in 
eq.(2.2) tends to sin (f) 2 (s') s'^ 1 /(s' — s) for s' — > +00, with \v 2 \ << 1, and then the dispersion 
integral converges. For / = we have that |Oo(s) | behaves as s~ Ui)+1 / s\ for s — > +00. We know 
already (f>o(s) for s quite high, and noticing that 0o(sff)/ 7r — 2 we take this value as a good 
approximation for v so that Qo(s) oc and o (+oo) = 2n. The integrand in eq.(2.4) tends then 
to sin0o(s')/( s/ — s) and converges because sin0 o (oo) — > 0. As a soft extrapolation for 4>o(s) for 
s > sh we shall take the same expression as in ref.[30] with 2n as the limiting value, 

MS) = ^ (l ± ,8>8 H , (2.13) 

with A 2 e [0,0.35] GeV 2 . In this way at s = sh one has an uncertainty of ~ ±7r/2 degrees that 
narrows only logarithmically for higher energies. This uncertainty is incorporated in the error 
analyses that follow. We think that our choices for (f>i(s) for s > sh, and the uncertainties taken, 
are quite reasonable. Let us note that for s < 1 GeV 2 the integrand in the dispersive integral is 
overweightly dominated for s' « sh- 

3 Dynamical approach to 77 — > tttt. 

In this section we review and update the approach of ref.[19] based on UxPT. This approach 
extends the one loop calculation of 77 — > tt°tt°, 7r + 7r~ and K + K~ , implementing unitarity 
in coupled channels, in order to account of the meson-meson final state interactions, and by 
including vector and axial-vector resonance exchanges. For the derivation of the unitarization 
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formulae that we use below we refer to the original ref. [19], see also refs.[22,23] for later advances 
in the formalism. In the dispersive method discussed in the previous section this is accomplished 
by the use of the Omnes functions Qo(s) and cj 2 (s) in the dispersion relations eqs.(2.4) and (2.2), 
respectively. 

On the other hand, we also offer here the expressions for the Born terms and the crossed 
exchanges of vector (1 ) and axial- vector (1 ++ , l + ~) resonances, employed both in the dispersive 
and dynamical approaches. The exchange of vector resonances calculated from chiral Lagrangians 
was performed in ref. [64], while the inclusion of the 1 ++ axial ones was undertaken in ref. [65]. 
Here we have calculated these tree level graphs but our results are different to those of ref. [65] 
for the latter diagrams. Reassuringly, we agree with the expressions given in ref. [14] for the 
l + ~ axial-vector resonance exchanges, which are also included here, and that share the same 
Lorentz structure as the 1 ++ resonance exchange contributions. It is worth mentioning that the 
1 ++ resonance exchanges contribute already at 0(p 4 ) in x?T for 77 — > it + it~, while the other 
resonance exchanges start one order higher in 77 — > tttt. The formulae corresponding to the Born 
term and the different resonance exchanges are given below. 

3.1 Formulae for Fn and Fc 

The expressions for F N (s) and F c (s) of ref. [19] are enlarged now to include the exchanges of the 
K* and the l + ~ axial-vector nonet: 

S S S ( — - — ■ \ 

f n 

(^"X 71 " ^p7r+7r _ tAw+w~ tcir+n~^) tw+w^ ,7r°7r° 
+ (t x K + 2tK*K+K- + tAK+K- + ^CK+K-) ^K+K- ,7r°7r° , 
Fc = T B + T p ^ +7r - + t S A ^ +qT - + t c ^ +1T - + [tp^O^O + t^O^O + tcTrOyrc) ^0^0^+^- 

+ (t X K + 2t K * K + K - + t AK + K - + 5t CK + K -^j t K + K - t7T + 7T - . (3-14) 

The different terms in this equation are given below. Although not explicitly indicated in the 
subscripts, the contributions from amplitudes with the intermediate K°K° state are encoded 
in the brackets proportional to tx+K-^n, because this S-wave strong amplitude is the same as 
tK K ,u-K an d kaon masses are all taken the same (isospin limit). 

In the previous equations the strong S-wave amplitudes can be obtained from those with isospin 
defined as, 







2 ,/=o 

^ 7T7T,7T7V 


+ 


2 7r7r,7T7r 


^7r°7r ,7r°7r 




2 ,/=o 

^ 7r7r,7T7r 


+ 


2 7T7r,7r7r 


t-ir+iT- ,7r°7r° 




2 ,/=o 

^ 7r7r,7T7r 




_j.I=2 

^ 7r7r,7T7r 



tK+K-,ir°n° ~ t K + K - i7V + n ~7^KR,tttt • (3.15) 
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In ref.[19] the 1 = and 1 = 2 S-wave amplitudes were calculated as in ref. [1]. For I = 2 we 
employ directly the experimental phase shifts t^ nw = — 5^/f e «<5 w (s) 2 s i a S n (s) 2 . 

In eq.(3.14) the amplitudes T^y^o, T P)Vr7r , T K * K x are the crossed exchange of the u, p and K* 
vector resonances, in order, denoted generically by T VP p. On the other hand, T APP and T CPP 
refer to the exchange of the 1 ++ and l + ~ multiplets of axial-vector resonances, respectively. In 
addition, t xP represents the one loop %PT function, 



t xP (s) 



2e 2 
16tt 2 



1 + 



ml 



, 1 + <T P (s) . 
1 - cr P (s) 



(3.16) 



with m P the mass of the pseudoscalar P, either n or K, and ap(s) = \J\ — Amp/ s. Finally, t RP p 
corresponds to the unitarity loop with the pseudoscalars P, P and the crossed exchange of the 
resonance R, either vector or axial, calculated as in ref. [19]. 

Here we give the expressions calculated for the Born terms, fig.l, 7(^1, ^i) / y(k 2 , A 2 ) — > P + (pi) P~ (p 2 ) , 



2e 2 



ei£ 2 - 



(eiPi)(e 2 p 2 ) (e 1 p 2 )(e 2 Pi 



Piki — i0 + P\k 2 — i0~ 



(3.17) 



where e\(ki, Ai), e 2 (k 2 , A 2 ) are the polarization vectors of the photons taken in the gauge ei&4 = 
ei/c 2 = e 2 /ci = e 2 k 2 = 0, with polarizations given by Ai and A 2 , respectively. In the following, we 
shall only need the three- momenta ki and k 2 along the z-axis and use the four- vectors (0, e(k, A)) 
with e(kz, ±) = T^(x ± iy) an d e[—kz, ±) = T^ijj ± ^r). The projection of the Born term in 
S-wave is given by i f dpT p {\i = +, A 2 = +,p)/4ir, with p the CM three- momentum of the pions. 
The prefactor % is included for convenience. An analogous formula for the projection in S-wave is 
employed for any other amplitude of 77 — > PP. We then have, 



Tf = 2e 2 ^log 1 + ^ (s) 



a n (s) 1 - a n (s) ' 



(3.18) 



with a w = \J\ — Am\j s. 

We now consider the crossed exchange of the vector resonances in 7(^1, Ai)7(/c 2 , A 2 ) — > P(pi)P(p 2 ). 
We evaluate them employing the chiral Lagrangian [64], 



C VP , = Ru,e e^F^Tr [$ {Q, d a V?}} , 



(3.19) 



where the matrix fields 3>(x), V(x) contain the octet of pseudoscalars and vectors [64] and Q is 
the quark charge matrix, Q = diag(2/3, —1/3, —1/3). 



V,PP 



e 2 R 2 v 



(eie 2 ) [s(t + nip) - (t - m 2 n ) (u - m\ )] - 2s(eij?i) (e 2 pi) 
t - M v + i0+ 

(eie 2 )[s(u + nip) - (t - ml)(u - ml)] - 2s(e 1 p 1 )(e 2 p 1 ) 



(3.20) 



If — My + ^0 + 

with P = vr°, 7T+, if+ or if°and V = u, p or if*. Furthermore, i^ = R\» = ij£/9. In these 
expressions M u , M p and M^* are the masses of the a;, p and if*, respectively. The coupling R^ 
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is determined from the width T(u — > 77r°), 



a (Af2 - ml) 3 
The projection of eq.(3.20) is given by, 



6 M 3 

= - 7775 " 2 v* r (" - 7vr°) = 1.44 GeV- 2 . (3.21) 



Typp — 2e Ry 



My . 1 + Op + Sy/s 

— log i — r ~ s 

Op 1 - (Jp + Sy/ S 



(3.22) 



with sy = 2{My — mp). It is important to remark that eqs.(3.20) and (3.22) contribute at the 
level of a two loop %PT calculation in 77 — > 7r7r. This implies that vector exchange is relatively 
suppressed at low energies due to its rather high chiral counting. Nonetheless in fig. 4 it is shown 
that they are significant above 0.4 GeV. The 1 ++ axial-vector crossed exchanges, as we shall 
see immediately, appear already at the level of one loop x?T. Hence, for low energy 77 — > Tin 
their contributions are more important than those due to the exchange of vectors. To evaluate 
the contribution from the crossed exchanges of the 1 ++ axial resonances in 77 — > 7r + 7r~ and 
77 — > K + K~ we consider the Lagrangian [16], 

cFa 

£ap, = i-^F^Tr ($[Q, VI) » (3-23) 

with a matrix with the 1 ++ octet of axial resonances and f n the pion decay constant. Here the 
antisymmetric tensor formalism for the axial-vector field is employed. The resulting amplitudes 
are 

rp p 2 F\ f, , l (e l e 2 )(t-mlf-2s(e lPl )(e 2Pl ) 
A '^~ M\fl I ( 1 2) " 4 M%-t-iO+ , , 

4 M 2 A -u-iQ+ J ' 

and the same expression for T AK + K - replacing — > m K . We use Ma = 1-23 GeV ~ \f2M p 
because of Weinberg's sum rules [66] and the KSFR relation [67]. We do not coincide with the 
expression given in ref.[65] (also employed in ref. [18]) for Ta^+tt-- However, we do agree with the 
formula given in ref. [14] for the similar exchange of the l + ~ axial- vector resonances, that has the 
same Lorentz structure. In ref. [65] the authors perform the check that at 0( P 4 ) they reproduce 
the xPT result for 77 — > 7r + 7r~ . This implies the axial- vector saturation of the chiral counterterm 
Lg + 1/ 10 , 

L 9 + L10 = t^j • (3-25) 

F 2 

We also fulfill this condition since at O ( P 4 ) our expression is Ta tt+tt- = 2e 2 M s A f2 (kik 2 )(e\e 2 ) + 0( P 6 ) 
and agrees with that of ref. [11] if eq.(3.25) is taking into account. We use eq.(3.25) to fix Fa from 
the value L 9 + L 10 = (1.4 ± 0.3) ■ 10 -3 . In the following we shall replace F\ by L 9 + L 10 using 
eq.(3.25). In ref. [68] the couplings of the vector and axial- vector resonances to pseudoscalars and 
external sources are calculated within the Extended Nambu-Jona-Laisinio model [69]. 
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The S-wave projection of eq.(3.24) is, 

2 L 9 + L w 



T 



A,PP 



2e A 



M\ 1 + a P 
log 



s A /s 



f 2 [<Jp 1 - Op + S A / s 



(3.26) 



with P = 7T+, K + and s A = 2(Mj 



m. 



). This equation differs with the S-wave projection of 



ref. [18], that employs the results of ref. [65] which disagree with our eq.(3.24). While we have the 
coefficient M\ in front of the logarithm, ref. [18] has s^/2. This mistake was also performed in 
ref. [19], because the expression from ref. [18] was taken. In addition, we have also included now the 
exchange of the K* vector resonance, not taken into account in ref. [19]. The numerical differences 
of correcting for these latter deficiencies are however very small. 

We finally give the expressions for the exchange of the l + ~ nonet of axial-vector resonances. 
We make use of the lowest order Lagrangian, see e.g. ref. [14], 



eD 

£=—F, u (C»{Q,d^}) , 

Jit 

where = -^C^X 1 + -j^C® ■ 1, with standard notation. We then have, 



T, 



C,PP 



with £ w o 



e 2 D 2 
18/?' 
= 5, ^ 



(eie 2 )(t 



m 



2 \2 



Ml 



2s(e ipi)(e 2 pi) (eie 2 )(tt 



m 



2 \2 



2s(eipi)(e 2 pi) 



u 



Ml 



1/2, £ K + = 1/2 and £ K o = 2. The S-wave projection is 



T, 



s 

C,PP 



e 2 D 2 



s + 



M, 



a 



log 



1 - Up + s c /s 



(3.27) 



(3.28) 



(3.29) 



Op 1 + Op + s c /s_ 

We take for Mq = 1.23 GeV, the value corresponding to the mass of the &i(1235) resonance. The 
constant D is determined from the bi — > jn + decay [70] , 

Mlr^ - 77T+) 216/2 



D 



(M 2 C - ml) 



~ 3.2-10" 



a 



(3.30) 



3.2 Isospin amplitudes 

Here we present our isospin amplitudes for 77 — > We employ the isospin states, 

-1 



F7T/0 



H + H + 7T 7T + + TV TV ) , 



0^0\ 



|7T7r) 2 = —=\7l + 7l +71 7T + — 27T7T' 

V6 



rpS 

J 77,0 



rpS 

J 77 ,2 



2 (rpS , X T 5 



(3.31) 



with the subscript indicating the corresponding isospin. As a result, for the S-wave of 77 — > (7r7r)/ 



(3.32) 



This isospin decomposition applies not only to the whole S-wave amplitude but to all of its partial 
contributions, Born term, vector and axial exchanges, etc. 
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4 Results 



In this section we discuss about the results that follow from the use of the dispersive and dynamical 
approaches, sections 2 and 3, respectively. 

4.1 Determination of the constants in the dispersive method 

Let us first consider the dispersive approach. The Li(s) functions have contributions from the 
Born term and from the vector and axial-vector crossed exchanges: 



2 / „e „c „c 1 



0) \ Tb + Tf j7r7r — Tj i7r7r + -^Tq^k ) 5 (4.33) 



where 



e^D 2 

"72" 

J 7T 



1 - (T V + S C /S 



s + — log 



(4.34) 



(Jjr 1 + OV + Sa/s. 

We now turn to the coefficients Cq and C2 of eq.(2.4). As discussed in section 2 these coefficients 
are determined by matching with the xPT one loop results for the linear behaviour in s of 
Fc( s ) — Bc(s) and Fn(s) when s — > 0. In order to impose the condition 1 one should realize that 
only T| 7r+?r _ is 0(s), being all the other resonance exchanges 0(s 2 ). Then, 

F c (s) = T§ + ^x^^s - r f 3 -f 6 + 0(s>) , (4.35) 
where x = (1 - m 2 /2M^)/(l - ml/M 2 A ). The one loop %PT result is 

and it follows that 



F c (s) = Ti + 4e 2 ^f^s + 0(s 2 ) , (4.36) 



with 



c + i=c 2 = f (4.37) 

e = 4v / 3e 2 ^±^( X -l) . (4.38) 

Notice that x is ver Y close to one and then £ ~ 

Regarding the condition 2, the linear dependence in s of Fjy(s) around s = from eq.(2.4) is 
given by 

^(s) = --^c s + y|c 2 s + 0(s 2 ) . (4.39) 
One loop x?T establishes that 

= -e 2 s/9Qn 2 f w + 0(s 2 ) , (4.40) 
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so that 



e 2 y/I73 2 

We evaluate the previous expressions given by one loop either by employing f n = 92.4 MeV 
or / ~ 0.94/^, where the former is the pion decay constant and the latter is the same but in the 
SU{2) chiral limit [7]. This amounts around a 12% of uncertainty in the evaluation of c and 
c 2 , due to the square dependence on f w . Note that both choices, f w or /, are consistent with 
the precision of the one loop calculation and the variation in the results is an estimate for higher 
orders corrections. This is included in the error analysis. 

The constant F (si) — L (si) in eq.(2.4) is related with the position of the Adler zero sa, where 
F N (s A ) = 0, by 



[F (si) - L (si)] = V2L 2 (sa) - L (sa) + (v^^sa) - c Qo(s A )^ s A 
fKs\r, , x f°° L 2 (s') sin <f) 2 (s')ds' 

4 r Lq(s') smMs')ds' 

« Wa) J«* s<\s>-s a )W)\ ■ 

As discussed in the condition 3 of sec. 2 the value of F (si) — ^o( s i) is finally fixed by requiring 
that (7(77 — > 7r°7r°) < 40 nb at the /o(980) peak, at ~ Si. This is based on the experimental 
data for 77 — > 7r°7r° of ref.[13], see fig. 7 where the data including the /o(980) are shown, whose 
maximum value at the / (980) peak of (7(77 — ■> 7r°7r°) < 40 nb, errors included. In addition, we 
also take into account the resulting width r(/ (980) — > 77) = 205t^(stat)tnl(sys) eV from the 
Belle Collaboration [31], and then impose that the width calculated from our approach lies within 
the interval 60-380 eV, as follows from the previous value. This interval is compatible with the one 
recently given for the width r(/ (980) — > 77) in ref.[71] and with the value r(/o(980) — > 77) = 
0.2 KeV predicted by the dynamical approach [19]. Note that r(/ (980) — > 77) is proportional to 
the maximum value of (7(77 — > tt°-k°) at the / (980) peak. For practical purposes both constraints 
are equivalent. 

4.2 The 77 — > 7r°7r° cross section 

In our normalization, the cross section in terms of F N (s) is given by 

a( 7 7^vrV) = ^-|^( S )| 2 . (4.43) 

Since the main contribution in the low energy region to f2/(s) comes from the low energy nn 
phase shifts one needs to be rather precise for the input parameterizations of the irn scattering 
data in this energy region in order to apply the dispersive method, eqs.(2.2) and (2.4). The small 
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V? (MeV) 



Figure 4: The 77 — > 7r°7r° cross section depending on the contributions included for Lj(s). Dot-dashed 
line (Born term only), dashed line (Born term plus vector crossed exchanges) and dotted line (Born term 
plus axial- vector crossed exchange mechanisms). The solid line is the full result. 

1 = 2 S-wave mi phase shifts, which induce small final state interaction corrections anyhow, can 
be parameterized in simple terms and our fit compared to data can be seen in fig. 2. For 1 = 
one of the parameterizations that we consider is the K-matrix of Hyams et. al. [52], 

Kij(s) = aiOij/ix! - s) + (3ipj/(x 2 - s) + ^ , (4.44) 

where 

x{ /2 = 0.11 ±0.15 
c*i = 2.28 ±0.08 
A = -1.00 ±0.03 
7n = 2.86 ±0.15 

in units of appropriate powers of GeV 

T = Sny/siK- 1 -iQY 1 , (4.46) 

where Q = diag(q. K ,qK) with q^qx) the center of mass momentum of the pions (kaons). The 
central values in eq.(4.45) correspond to S n (sK)o < 7r. However, varying the parameters within 
their errors it is possible to obtain ^(5^)0 > ^ (see fig. 2 where both cases are shown). Since 
below 0.6 GeV this K-matrix was not fitted to data and strongly disagrees with K e 4 data [55, 56], 
we consider instead the parameterizations of ref. [59] (CGL), ref.[60] (PY) and ref. [10] (AO). They 
agree with data from K ei decays [55,56] and span to a large extent the band of theoretical 
uncertainties in the 1 = S-wave irn phase shifts [52-54]. PY runs through the higher values 
of S n (s) , while CGL does through lower values, see fig. 2. AO, double-dash-dotted line, moves 
somewhat in between the previous results for energies between 0.5-0.7 GeV and then follows closely 



xl /2 = 1.19 ±0.01 

a 2 = 2.02 ±0.11 , v 

I3 2 = 0.47 ±0.05 v ' ; 

712 = 1.85 ±0.18 7 22 = 1.00 ±0.53 , 

. The K matrix is related to the 1 = S-wave T— matrix by 
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CGL results. CGL is used up to 0.8 GeV, since this is the upper limit of its analysis, and the K- 
matrix eq.(4.44) above that energy. On the other hand, PY is employed up to 0.9 GeV, as in this 
energy this parameterization agrees well inside errors with [52], and eq.(4.44) above it. AO is used 
up to 1.5 GeV since it fits S-wave meson- meson data up to 2 GeV. We shall also employ in some 
cases the results for the tttt phase shifts from ref. [1], which gives rise to a /o(980) signal tested in 
many other processes like e.g. J/\I/ [50,72], D [73] and [20,74] decays. This parameterization 
will be referred in the following as 00 and will be used in in the region 0.9 < y/s < 1.1 GeV, 
above it eq.(4.44) is again taken and below CGL. 

140 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 'J 





Vs (MeV) 



Figure 5: I = (dashed) and 2 (dot-dashed) contributions to the 77 — > 7r°7r° cross section. The solid 
line is the full result. The curves have been generated employing the parameterizations CGL and Hyams 
[52]. 

In fig. 4 we show the different contributions to the 77 — > 7r°7r° cross section compared with 
data* 3 [13]. The different curves originate by considering only the Born term as the only contri- 
bution to Lj(s) in eq.(2.4) (dot-dashed line), the Born term and the exchanges of vector resonances 
(dashed line), the Born term and the exchanges of axial- vector resonances (dotted line) and the 
full model (solid line) for L Iy i.e., the Born term plus vector and axial- vector resonance exchanges. 
The curves have been evaluated with the central values of the parameterizations eq.(4.44) and 
CGL. One can see that the dominant contribution comes from the dressing of the Born term by 
final state interactions, despite that the Born term at tree level is zero. The vector and axial-vector 
resonance exchange mechanisms also give important contributions to the total cross section, even 
already at threshold for the exchange of the axial-vector resonances. Most of the vector resonance 
crossed exchange contribution comes from the u exchange, since the u coupling to 771"° is about 
one order of magnitude larger than that due to the p or K* resonances. It is clear from the 
figure that the axial- vector meson exchange mechanisms, overlooked in refs.[26, 28], are found to 
be substantial and tend to increase the cross section. Notice that the cross section given in ref. [26] 
is very close to the dashed line in fig.4. 

# 3 Since the data only cover up to | cos#| < Z, with Z — 0.8, we have divided them by the factor Z. This suited 
for low energies as S-wave dominates, for higher energies is a normalization convention. 
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It is also worth stressing that there are strong interferences between the different mechanisms, 
as well as between the isospin 1 = and 1 = 2 amplitudes. In fig. 5, we show the 77 — > 7r°7r° 
cross section considering only the 1 = S-wave contribution (dashed line), the 1 = 2 S-wave 
contribution (dot-dashed line) and the full result (solid line). For definiteness, the curves are 
calculated with the central value of ref.[52] and CGL for low energies. We can see the strong 
cancellation between the 1 = and 1 = 2 contributions due to the absence of the Born term at 
tree level in the 77 — ► 7r°7r° channel, which is the dominant contribution in the 1 = and 1 = 2 
channels. The Born term is the responsible in every isospin channel for the bumps at low energies 
in fig. 5, though the strong final state interactions in I = reshapes it considerably and tends to 
accumulate strength to lower energies than for the 1 = 2 part. 




Figure 6: The error band around the solid line is the uncertainty in our results by requiring that 
17(77 —* 7r°7r°) < 40 nb at s\. The solid line should be compared with the dot-dashed one, that would 
correspond to the formalism of ref.[26] in the case that (fio(s) ~ S n (s)o above sk and including the 
axial- vector exchanges. The lines are evaluated with the central values of CGL and eq.(4.44). 

In fig. 6 we show with the solid line our results for the central values of eq.(4.44) and CGL while 
the dashed-solid line corresponds to the result that would be obtained if using the approach of 
ref. [26] (including also the axial- vector exchanges). The improvement achieved by passing from the 
dot-dashed line to the solid one is determined by using eq.(2.4) instead of eq.(2.2), being the latter 
the original equation derived in refs. [26-28] . The gray band around the solid line corresponds to 
the uncertainty originated exclusively by taking 17(77 — > 7r°7r°) < 40 nb at s\. 

We include in fig. 7 the /o(980) region as well, plotting for y/s < 1.05 GeV. The generation 
of the /o(980) resonance peak within the dispersive method is a remarkable achievement with 
important consequences as it has been used already for fixing that the z = +1 behaviour of 4>o(s) 
is the one realized in nature. We have also used it to constraint the values of F (si) — L (si), 
so that we are able to take three subtractions in the dispersion relation instead of just two, as 
in previous studies. The values of F (si) — L (si) are fixed in this figure such that the size of 
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Figure 7: The 17(77 ~~ > ^ ^ ) cross section including the region of the /o(980). We show two different 
curves such that the cross section at the top of the peak is 10 and 30 nb, long-dashed and solid lines, 
respectively. The short-dashed line employs the 00 parameterization around the /o(980) resonance. The 
dot-dashed line is evaluated with the dynamical method of ref.[19]. 

<t(77 — > 7r°7r°) in the / (980) peak is 10 or 30 nb, dashed and solid lines, respectively. These lines 
are evaluated with the central values of CGL and eq.(4.44) for definiteness. From this figure it 
is clear that above ~ 0.8 GeV the D-wave contribution in 17(77 ~~ * 7r°7r°), due to the tail of the 
/ 2 (1270), is important and this resonance is needed so as to agree width data. This contribution is 
added incoherently because of the azimuthal integration (the ^2(1270) resonance signal in 77 has 
dominantly helicity 2 [27]). We also show in the figure by the short-dashed line the result from 
the dispersive approach but employing the parameterization OO for the tttt phase shifts around 
the /o(980) resonance. The dot-dashed line is calculated with the dynamical method of ref. [19]. 
This method rightly predicts the peak at 0.5 GeV, with size very similar to that obtained from 
the dispersive method and using PY, though tends to give somewhat a too fast decreasing cross 
section for energies above > 0.6 GeV. In fig. 9 we also show the curve corresponding to the AO 
parameterization with a peak value in the /o(980) of 20 nb. In this way one can see how the actual 
shape of the /o(980) signal changes depending on the parameterizations and on the peak values 
used. Further theoretical precision in this energy region requires to include the D-wave 77 — > tctc 
contribution and fit simultaneously the recent data on 77 — > 7r + 7r~ [31]. This is beyond the scope 
of the present study. 

Our final 7(77 — > for y/s < 0.8 GeV is shown in fig. 8. We give the results for CGL(solid) 

and PY(dashed) and eq.(4.44), as usual. The corresponding band around the lines has into account 
all the theoretical uncertainties discussed with respect to the dispersive approach. For that, we 
implement a Monte-Carlo Gaussian sampling taking into account the errorbars in the parameters 
of CGL, PY and eq.(4.44), the uncertainty for s > sk in the approximation 4>o(s) ~ 5 7r (s) in 
the region 1 < y/s < 1.5 GeV and that from eq.(2.13) above sh- We also include in the error 
the 12% uncertainty in evaluating c and c 2 according to eq.(4.41) and the variation within the 



20 




son 



Figure 8: Final result for the 77 —* 7r°7r° cross section with the theoretical error band. The experimental 
data are from the Crystal Ball Collaboration ref.[13]. The dashed line is evaluated with PY and the solid 
line with CGL. The band around each line corresponds to the theoretical uncertainties. The dot-dashed 
line corresponds to AO and does not include error band, of similar size than the ones shown, in order 
not to overload the figure. 



bound a(77 — > 7r°7r°) < 40 nb at s\. The width of the band is dominated by the errors in the 
parameterizations PY and CGL. In the same figure the dark dot-dashed line corresponds to the 
AO phase shifts. We do not show its corresponding error bar as it is of similar size as for CGL 
and PY and it would complicate the visualization of the information contained in the figure. In 
addition, the dotted line corresponds to one loop CHPT [11,12] and the lighter dot-dashed one 
is the two loop result [14,15]. The latter is closer to the dispersive results for low energies but 
still one observes that the 0(p s ) corrections would be sizable. It is worth stressing that if the 
axial- vector exchanges were removed, in refs.[26, 28] they were not included, then our curves would 
be smaller. This corresponds to the dot-dot-dashed line in fig. 8, which is very close to that of 
ref. [26] when employing 0o( s ) — <5tt( s )o — 7r for s > sk- The former curve is evaluated making use 
of CGL and ref. [52] with their central values. The three experimental points [13] in the region 
0.45-0.6 GeV agree well with this curve. However, once the axial-vector are included the curve 
rises. These three points lie around 1.5<t below the CGL result band, and by somewhat more 
than two sigmas below the PY band. This clearly shows that more precise experimental data 
on 77 — > 7T°7T° could be used to distinguish between different low-energy parameterizations of the 
/ = 7T7T S-wave phase shifts. The next three experimental points in the region 0.6-0.75 GeV agree 
better when the axial-vector resonance contributions are taken into account, as one should do. 
As a result of this discussion, more precise experimental data for 77 — > 7r°7r° are called for. The 
width of the error bands around the solid and dashed lines is reduced as compared with ref. [29] 
because now the z = +1 case has been singled out. 

In fig. 9 we show our final results for 17(77 ~~ > 7r°7r°) including the /o(980) region as well. The 
error bands are evaluated similarly as for fig. 8. One appreciates how the widths of the error bands 
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Figure 9: Final result for the 77 — ► 7r°7r° cross for \fs < 1.05 GeV. The experimental data are from the 
Crystal Ball Collaboration ref.[13]. The dashed line is evaluated with PY and the solid line with CGL. 
The band around each line corresponds to the theoretical uncertainties. The dot-dashed line correspond 
to using AO. Its error band, not shown, is similar to the ones of the other two curves. 

largely increase in the /o(980) region because of the uncertainty in the input for the value at 
the peak 4 < (7(77 — > 7r°7r°) < 40 nb, which is obtained by requiring that r(/o(980) — > 77), as 
calculated from our approach, lies in the range of the value determined in ref.[31], r(/ (980) — > 
77) = 205ts3(stat)t{il(sys) eV. This bound in the size of the / (980) peak could be also inferred 
by simple inspection of the experimental data in the same figure since otherwise the /o(980) peak 
would be larger than the total cross-section, which also includes incoherently a large D-wave 
contribution from the tail of the /2(1270). A reduction of this uncertainty in the /o(980) has to 
await a proper study of the new data 77 — > -k + -k~ [31] along these lines, as already commented. 

From our approach we obtain that the Adler zero in F^{s) is located at sa = (1.11 ± 0.03) m\, 
if CGL is used, and at sa = (1.09 ± 0.04) m^., when PY is taken for lower energies. If AO is used 
up to 1.4 GeV we obtain for sa the same value as the latter one. From our analysis then sa lies 
in the range (1.05-1.14) m\. This range is in agreement with the the two loop value [14], 
sa = 1.175 m^. The difference is around a 4%, well inside the expected uncertainties of ref.[14]. 

4.3 The a — > 77 width 

Let us recall that the a is the lightest resonance with the quantum numbers of the vacuum. It was 
originally included by studies of the AW potential, but due to the advent of the chiral potentials 
[75-77] it is no longer necessary and, furthermore, chiral symmetry seems to substantially reduce 
its contributions [78]. Another relevant example where the a is introduced is the linear sigma 
model [79]. However, when the non-linear sigma models are unitarized, fulfilling the unitarity 
and analyticity requirements associated with the right hand cut, a pole associated with the a 
resonance appears as well [1-5,80]. The works [1,3] were the first to show that the presence of 
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such a light scalar resonance is not in contradiction with the non-linear sigma models. It appears 
because of the strong interactions, driven by chiral symmetry, between the two pions in an / = 
S-wave pair. In ref.[81] it was shown that in the limit of chiral symmetry restoration, f n — > 0, the 
a pole moves to zero and becomes the chiral partner of the pion. In this respect, the difference 
in mass between the a resonance and the it is an order parameter for chiral symmetry breaking. 
It is interesting to consider the a coupling to 77 because, as indicated in the Introduction, it is a 
complementary information apart from that of scattering as it is sensitive to the electric charge. 
This can give clues to possible underlying structures in its nature. 

For the evaluation of the coupling a — > 77, g all1 the amplitude F^{s) has to be evaluated not 
on the first or physical Riemann sheet but on the second one. This is accomplished by shifting 
q n — > —q n . We denote by F N (s) the resulting amplitude on the second Riemann sheet. The strong 
I = nn S-wave T— matrix on the same sheet, T n (s), can be calculated in terms of the one on 
the physical sheet, T/(s), by 

T n and F N have a pole corresponding to the a resonance at s a , so that 

lim Fx = ^ 77WV ° , (4.48) 

5 >Sa Sfj S 

where g an o n o is the coupling of the a to 7r°7r°. #4 

From F N of eq.(3.14), corresponding to the extended dynamical approach of ref.[19], the cou- 
pling g a ~ n is 

5 , ct77 — Y "(^Avr+vr- + l-l^Cir°7r + ^W7r°ir° + 2t p7r + 7r - + t X7T )g an7T 



+ -^J={pAK+K- + 5tcK+K- + 2tK*K+K- + ^xK)9ctKK 

Y^^Att+tt- + l-l*Cir°ir° + *a;ir 7r° + 2t pn + n - + t 
H 7= {tAK+K- + 5tcK+K- + 2t K *K+K- + t x x) ^° 



(4.49) 



V2 go 

In the previous equation g U7Tn is the coupling of the a to the I = mr state in the unitary 
normalization of refs.[19, 80]* 5 and, analogously, g a Kk is the coupling to the 1 = KK one. 
Notice that the ratio of couplings is given directly by g^Kkl AWtt — t° Kn > kkI^-k-k kk' w ith the 
strong amplitudes evaluated on the second Riemann sheet. 

It is important to remark that the ratio #0-77/ go-™ does not depend* 6 on the strong meson-meson 
coupling, hence it allows a clearer comparison between different approaches for the 77 — > mi part. 



# 4 As it is well known [82] the phase of a resonance coupling is background dependent. Hence, the phase of gan"^ 
extracted from 77 — > ir°ir° differs from that obtained, e.g., from 7r°7r° — > 7r°7r°. This does not bother us here since 
we are only interested in the moduli of g ail and j ffT o T o. 

# 5 With the unitarity normalization the width to tttt of an I = narrow resonance of mass M is given by 
r = |<w|V/(8ttM 2 ). 

# 6 Actually there is a dependence from eq.(4.49) on the strong coupling due to the kaons, but this term is very 
small numerically. 
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We now determine g a ~ n /g aiT7T and by using the dispersive approach of section 2. The 

dispersion relation eq.(2.4) to calculate F (s) is only valid on the first Riemann sheet. If evaluated 
on the second Riemann sheet there would be an extra term due to the a pole. However, the relation 
between F and F , the latter on the second sheet, can be easily established using unitarity above 
the 7T7T threshold, 

F Q (s + ie) — F (s — ie) = —2iF (s + ie)p(s + ie)T 7 ° 7 (s - ie) , (4.50) 

with 4m\ < s < 4m|-, p(s) = a 7T (s)/167T and e — > + . Due to the continuity when changing from 
one sheet to the other, 

F (s - ie) = F (s + ie) , T/ =0 (s - ie) = T}f°(s + ie) . (4.51) 

it results then from eq.(4.50) 

F (s) = F (s) (1 + 2ip(s)T I I f{s)) , (4.52) 
where the analytical extrapolation of eqs.(4.50) and eq.(4.51) has been taken. Around the a pole, 

T i=o = __gL L _ ; p Q ( s j = ^ g*~ n g*™ ^ (4 53) 

s a s s a s 

and then from eq.(4.52) 




(4.54) 



This expression allows to evaluate again the ratio g all / g a -mv 

We consider in the subsequent s a = (M CT — iT (T /2) 2 either from the studies of ref. [10] or ref. [8]. 
The former provides a reproduction of the S-wave 1 = and K~n + — > K~ir + data from tttt 
threshold up to around 2 GeV, considering for the first time 13 coupled channels for the / = 

meson-meson S-wave. The interactions kernels are derived from chiral Lagrangians and are 
unitarized employing the N/D method [2]. Multiparticle states are mimicked by aa, pp, ww, 
The aa elementary transition amplitudes are given without any new free parameter, the same as 
for the vector- vector ones which are derived by making use of a massive Yang-Mills theory [83]. 
In ref. [10] one has = (456 ± 6) MeV and = (482 ± 20) MeV. Ref. [8] provides the values 
m ccl = 441+16 Me y and t ccl = 544+is MeV. The latter values result from the / = tttt S-wave 
amplitude obtained after solving the Roy equations and matching with CHPT to two loops. In 
the following the superscripts AO and CCL refer to those results obtained by employing s a from 
ref. [10] or [8], respectively. 

In table 1 we give | g all / g a7rn \ for the dispersive and dynamical approaches of sections 2 and 
3, respectively. The second and third rows correspond to the dispersive approach with s a from 
ref. [8] and [10], respectively. The result from the dynamical approach are shown in the last row. 
This approach makes use of strong amplitudes with s a very close to that of ref. [80] . Notice that in 
order to calculate -F/(s), and then 17(77 — > 7r°7r°), s a is not needed. The different columns in table 

1 represent different contributions to Lj. From left to right in table 1: only the Born term is taken 
as left hand cut, Born term and vector resonance exchanges, Born term and axial- vector resonance 
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Qa-yy/ Qu-ktt X 10 


Born 


Born + vector 


Born + axial 


All 


Dispersive approach 
s a from ref . [8] 


2.60 


2.22 


2.37 


2.01 ±0.11 


Dispersive approach 
s fj from ref. [10] 


2.42 


2.02 


2.25 


1.85 ±0.09 


Dynamical approach 


2.05 


1.84 


1.81 


1.61 



Table 1: Different contributions to the ratio \ga-yy / g^nnl given by the dispersive approach of section 
2, second and third row, and the dynamical one of section 3, fourth row. In the second row the s a 
pole position is taken from ref. [8], while for the third one s a comes from ref. [10]. No error is given 
for the dynamical approach since this model is very constraint and the systematic uncertainties 
cannot be worked out. 



exchanges and all together. We observe that the interferences between the different pieces has 
opposite sign to those in the 77 — > it it cross section, fig.4. This difference occurs because in the 
7r°7r° cross section both the 1 = and 1 = 2 play a role, while in the a width only the 1 = takes 
part. For the dispersive results we have taken the average the resulting ratios \g al ^/ g a7Tn \ obtained 
with PY and CGL parameterizations. The values obtained if using AO agree with the values 
shown within one sigma. No errorbar is given for the dynamical model result because this model 
is very constraint, has no free parameters, and the systematic uncertainty cannot be estimated. 
Its value is quite close to that of the second column, with s a from [10], only a 10% lower. 

These ratios of residua at the a pole position are the well defined predictions that follow from 
our Fq(s), given s a . However, the radiative width to 77 for a wide resonance like the a, though 
more intuitive, has experimental determinations that are parameterization dependent. This is 
due to the non-trivial interplay between background and the broad resonant signal as stressed in 
ref. [29]. An unambiguous definition is then required and we employ the standard narrow resonance 
width formula in terms of g CT77 as determined from the residue at s a , 

as done in refs.[26, 29]. Nevertheless, the determinations of the radiative widths from this expres- 
sion and those from common experimental analyses of physical data can differ at the level of a 
< 20%, as discussed in ref. [29]. Notice that in the equations usually employed in phenomenological 
fits to data, e.g. see ref. [31], the couplings are determined along the real axis, while the definition 
eq.(4.55) requires to move to the pole position. 

The use of a narrow width resonance formula eq.(4.55) in order to evaluate T(a — * 77) is a 
precise definition and convenient criterion, though somewhat arbitrary as the a meson has a large 
width. Elaborating a bit more on this issue let us introduce the normalized a mass distribution 
D(s), 

/■oo 

D(s)ds = 1 . (4.56) 



[•CO 

Jiml 
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Then, 

! 2 /■ ^ 



|5 , o'77 



a— >77 



7JL> W - (4 ' 57) 



8nM a J Ami s 1 / 2 ' 

with g 7 = is the photon center of mass three-momentum. We can then rewrite eq.(4.57) as, 

^■■^kiC^-ti- < 4 - 58 > 

In the last step we have taken into account eq.(4.56). Despite the introduction of the mass 
distribution for the a, D(s), eq.(4.58) corresponds to the narrow resonance width approximation 
due to the null mass of the photons. Of course, in the previous equations we have assumed 
an energy independent coupling g ail . A more physical calculation would imply to take this 
dependence into account but this is beyond the scope of the present paper. In the following, we 
shall content ourselves with eq.(4.55) as a well defined translation from |<7o- 77 | to the more familiar 
concept of width. 

In order to determine T(a — > 77) by applying eq.(4.55) with the values given in table 1 one 
needs to provide numbers for | Po-tttt I - This is a crucial input to evaluate the width to two photons 
as the latter is proportional to \g a7T7r \ 2 . We first consider the value \g^ w \ = (3.17±0.10) GeV from 
the approach of ref. [10]. # 7 The calculated width is then 

T AO (a^nn) = (1.50 ± 0.18) KeV . (4.59) 

Not only the position of the pole in the partial wave amplitude, but also its residue can be 
calculated in the framework of the dispersive analysis described in ref. [8]. Expressed in terms of 
the complex coefficient g an7T defined in eq.(4.53), the preliminary result for the residue amounts 
to = (3.311°;$ GeV.# 8 

V ccl {(T ^ 77) = (1.98+o.so) KeV ? (4 60) 

Taking the average between the values in eqs.(4.59) and (4.60) we end with, 

r(o- -> 77) = (1.68 ±0.15) KeV . (4.61) 

This number is in agreement with the values (1.8 ± 0.4) KeV, with s a from ref. [80], and (2.1 ± 
0.3) KeV, with s a from ref. [8], calculated in ref. [29]. Note that we use here the most recent analysis 
of ref. [10] instead of ref. [80] for the a pole position and coupling. In addition, the coupling \g„^\ 
has been also explicitly evaluated from ref. [8], which was not known at the time of ref s. [26, 29]. 
On the other hand, eq.(4.61) agrees at the level of one sigma with the estimate T(a — > 77) = 
(1.2 ± 0.4) KeV of ref. [84], that considers the experimental data on the proton electromagnetic 
polarizabilities. 

Pennington in ref. [26] predicted T CCL (a — > 77) = (4.09 ± 0.29) KeV for s a from ref. [8]. 
The difference between this number and our result in eq.(4.60) is due to i) ref. [26] calculates 
\g ajl /g a7r7V \ = (2.53 ± 0.09) x 10~ 3 for s a from ref. [8] instead of (2.01 ± 0.11) x 10~ 3 as given in 



#7 We thank M. Albaladejo for supplying this number before publication. 

# 8 We express our gratitude to I. Caprini and H. Leutwyler for sending us this result for the residue. 
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the second row of table 1. This is due to the omission of the axial- vector resonance exchanges 
in ref. [26,28], 10%, and to the improvement in our approach by taking into account one extra 
subtraction and slight different input for the phases, another 10%. One then has the suppression 
factor (2.01/2.53) 2 = 0.63 for the width T(cr — > 77). ii) There is also the difference with regards 
the value of the \g^^\ employed. Ref.[26] uses \g^^\ = 3.86 GeV instead of 3.3 GeV, with the 
last number very close to our previous estimate of \g^ n L \ in ref. [29] (see eq.(3.14) of this reference). 
The additional suppression factor (3.3/3.86) 2 = 0.73 then raises, hi) Finally, ref.[26] includes and 
extra factor |/?(s CT )| = 0.956 as compared with eq.(4.55). Taking all these several factors into 
account one has 4.09 KeV -> 4.09 x 0.63 x 0.73 x 1.05 = 1.98 KeV, as in eq.(4.60). 

As it is stressed in ref. [26], the a resonance is far from dominating the 77 — > 7r°7r° reaction. 
Were this the case, then 

F N (s) ~ 9 -^l^l , (4.62) 
s a s 

from where 

Given the value of (1(77 — > tt°-k°) ~ 16 nb in the peak, one then concludes that the T(a — > 
77) ~ 0.39 KeV, which is about a factor 4 smaller than our calculation eq.(4.61), and one order 
of magnitude smaller than the value of ref. [26]. In this latter reference it was argued that the 
reason for this difference is the very strong destructive interference between the 1 = 2 and 1 = 
amplitudes, see fig. 5. However, a quantitative analysis, like that offered in ref. [26], is affected by 
the fact that the ratio \F /F N \ 2 depends strongly on the precise energy where it is evaluated in 
the energy region where 77 — > tt°-k° peaks. This is clearly seen in fig. 5. At 400 MeV the ratio 
\Fq/Fn\ 2 ~ 12 (one has to multiply by 3 the ratio of the cross sections of / = and 7r°7r° to 
correct for the Clebsch-Gordan coefficient in eq.(3.31)), while at 500 MeV is ~ 3/2. This figure 
also shows that peak of the 77 — > 7r°7r° cross section at around 500 MeV is not directly related to 
the a resonance, but a a result of the strong interference between the 1 = 2 and 1 = components. 
Note also that the latter lacks of a resonance structure around the a due to the Adler zero in 
-^/v(s)-* 9 The a resonance, present in the 1 = component, distorts the form of the 1 = Born 
term accumulating strength towards the threshold so that above ^Js ~ 0.5 GeV |-fo(s)| is much 
smaller than the 1 = 2 counterpart. 



5 Conclusions 

We have performed a theoretical study of the reaction 77 — > 7r°7r° making use of improved versions 
of the approaches of refs.[29] and [19]. The former is based on a dispersion relation for 77 — > nn in 
S-wave and definite isospin / similar to that of refs. [26-28] , but including one more subtraction for 
1 = 0. This constant is fixed by taking into account an additional constraint based on a lax bound 
for the size of the /o(980) peak in 77 — > 7r°7r° or, equivalently, for the value of r(/ (980) — > 77) 
from ref. [31]. In this way, ref. [29] was able to drastically reduce the uncertainty in the results 
of ref. [26] on the function phase (f)o(s') of F (s) taken above the KK threshold. Here, we have 

#9 See refs. [73, 85] for a discussion on the role of the Adler zeroes and a resonant shapes in two pion production 
processes. 
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been able to reduce it further so that it is shown that 0o( s ) should follow quite closely S (s) 
above the KK threshold. Then, in the advent of better experimental data on (7(77 — > 7r°7r°), 
our present results allow a sharper comparison with data in order to distinguish between different 
parameterizations for the low energy tttt I = S-wave. 

We have also updated the approach of ref. [19] adding the crossed exchanges of the K* and the 
nonet of the lightest 1 ++ axial-vector resonances, not included in this reference. The axial-vector 
exchanges of the l + ~ and 1 ++ resonances were not either taken into account in ref. [26, 28], but 
have been included here. The explicit expressions and calculation details for the different elements 
both in the dispersive and dynamical approaches are also given. In particular, the expressions 
for the Born terms and the crossed exchanges of the 1 , l + ~ and 1 ++ vector and axial- vector 
resonances, respectively, are collected. 

The role of the a and / (980) resonances in 77 — > tttt has been also addressed. For the latter 
resonance it was argued by continuity arguments that it should come up as a small peak, as actually 
seen recently in ref. [31] for 77 — > 7r + 7r~. Importantly, the fact that it is indeed a peak requires 
that the phase of F (s) should rapidly increased around the KK threshold, the so called z — +1 
scenario. Regarding the a resonance, we have calculated the ratio of the residua \g aJ7 / g an7V \ = 
1.91±0.07 taking into account the values of the a pole position, s a , of refs.[8, 10]. One should stress 
that the calculation of (7(77 — ► rnr) is independent of the value taken for s a and |<? CT 7r7r|, though 
they are essential inputs for calculating the width T(a — > 77). With \g^ n \ = (3.17 ± 0.10) GeV 
for s a of ref. [10] one has T AO (a -> 77) = (1.50 ± 0.18) KeV and with \g^ n L \ = (3.3llg;££) GeV, 
a preliminary value corresponding to the s a of ref. [8], T CCL (cr — > 77) = (1. 981^24) KeV. The 
average between the previous two values reads T(a — > 77) = (1.68 ± 0.15) KeV, that we take as 
our final result for the width. 
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